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On , to the Higgs bosons are very small or forbidden. However, loop corrections may 

' generate the gauge-symmetry-breaking mass splitting within the multiplet, which 

does not vanish in the large M limit due to the threshold singularity. We calculate 
^ ' the electroweak contribution to the mass splitting for a heavy fermion multiplet, to 

■ the two-loop order. Numerically, two-loop electroweak contributions are typically 

O(MeV). 



Abstract 

New heavy particles in an S\]{2)i multiplet, sometimes introduced in extensions 
of the standard model, have highly degenerate tree-level mass M if their couplings 
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1 Introduction 



In some extensions of the standard model, there are new heavy particles which belong 
to an SU(2)£,xU(l) multiplet F and have no, or very small, mixing with other particles. 
The masses of these particles are almost degenerate to a value M by the gauge symmetry. 
Although the spontaneous breaking of the 811(2)/^ xU(l) symmetry by Higgs bosons may 
generate mass splitting 5M among them, the tree-level mass splitting generally behaves 
as 5M ~ m\r/M and becomes very small for M ^ mw- This is especially the case for 
a very heavy fermion multiplet where tree-level renormalizable couplings to the Higgs 
bosons FFH are forbidden by symmetry. Some of the examples are the almost pure 
winos or higgsinos [IlEllSllllElEllTllH] in special parameter regions of the minimal 
supersymmetric standard model [9], SU(2)£, triplet fermions in Type III seesaw model 
for neutrino masses [TOl [11], and also models [121 [IS] where vector-like heavy fermion 
multiplets are added to the standard model by hand. 

In such cases, it has been known [2], HI [6l El [HI [121 [131 [H] that the dominant part of the 
gauge-symmetry-breaking mass splitting within the multiplet F comes from the radiative 
correction. Although the form of the mass correction strongly depends on models, the 
contributions involving electroweak gauge bosons V = (7, Z, W), shown in Fig. 1(a) for the 
one-loop, are common in a wide class of extended models. Since gauge symmetry breaking 
in this diagram comes from the squared masses (m^, m|) in the loops, one naively expect 
the 0{a2m'^r/M) contribution to the mass splitting. However, due to the singularity of 
the diagram near the threshold, at = ~ (M-f-my)^, 0{a2mw) contribution to the 
mass splitting appears, which does not vanish in the M ^ mw limit: Roughly speaking, 
it is "nondecoupling" . This mass splitting is phenomenologically interesting, especially in 
the case where the neutral component /° of F, either fermion or boson, is stable or has 
very long lifetime, and may be a candidate for the cosmological dark matter. In such a 
case, the loop-generated mass splitting between charged components f^{Q 7^ 0) of F and 
/° is crucial for estimating the rates of the /"^ + ■ decays expected at colliders, 

and also for possible resonant annihilation /°/° — > f'^f"^ — »■ VV for indirect detection 
off 

To evaluate the mass splitting within F to the next-to-leading order, we need two- 
loop calculation of the mass correction for the members of F. In this paper, we perform 
such calculation for the loop corrections by the standard model particles, generated by 
the electroweak gauge interactions of F. For simplicity, we concentrate on the SI]{2)l- 
breaking and "nondecoupling" part of the mass correction, which should be relevant for 
the mass splitting in the M ^ mw case. 

2 One-loop mass correction 

Since the electroweak contributions to the mass correction should be determined by the 
SU(2)ixU(l) representation of F, we work in the framework of the Minimal Dark Matter 
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model [13] , which has been proposed as a minimahst approach to the dark matter problem, 
for the fermion case. In this case, Dirac or Majorana fermions in an SI]{2)l multiplet F 
with 811(2)1, isospin / and U(l) hypercharge Y (and having no SU(3) color) are added to 
the standard model. The lagrangian is 

C = CsM + cF[-f^^D,-M]F, (1) 

where c = 1(1/2) for Dirac(Majorana) fermions, respectively. Note that the mass cor- 
rections presented in this paper are common to both types of fermions. Z?^ denotes 
SU(2)iXU(l) gauge covariant derivative for F. Since F has no direct couplings to the 
Higgs boson, the members of F, (with charge Q = I3 + Y , I3 = —I, —I + 1, . . . , J) 
have a common mass M at the tree-level. We assume that M is sufficiently larger than 
the masses of standard model particles {W, Z, top quark t, Higgs boson h), typically 
M = O(TeV) which is cosmologically favored in the Minimal Dark Matter model [13]. 
We also use approximation that all other particles in the standard model are massless. 
The pole mass Mp of at the two-loop order is given in terms of the self energy of 

S(p) = S^^(p2)^ + SM(p'), (2) 



as 



M - Sm(m2; 



^ 1 + S^(M2) 

= M - [MSiJ^ (M^) + (M^)] - [MSi^^ (M^) + (M^)] 

+ [MJ:^^\m^) + si^)(M2)][S^)(M2) + 2M^t^^\M^) + 2Mt^l}{M^)] 
= M + 6M^^^ +6M^^\ (3) 

Here ^k]m ^k]m ^^e one-loop and two-loop parts, respectively. The dot in Eq. ([3]) 
denotes the derivative with respect to the external momentum squared. The absorptive 
part of the self energy is 0{g^) and need not be considered here. Loop integrals are 
regularized by the dimensional regularization {D = 4 — 2e) with the MS subtraction 
scheme. 

The form of the one- loop mass correction 5M(^) is well known [2l[6l[Hl[IIl[12l[l3l[Tl[T5]. 
Abbreviating the factor 0:2/(477), it is expressed as 

6MW = [Cp - li)Xl^^ + sl^ih + YfX^^^ + clih - AvYfx'-p, (4) 

where Cp = I{I + 1), cw = cos9w = mw/mz, sw = sinOw, tw = tan6'vF, and 



X^^ = M 



M 



'2 + ^ j Bo{M', M, my) - 1 + -^{A{M) - A{mv)} 



--31ogiVr + 4-/ 



M 
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>/2 — X 

f{x) = 2x{2 + x^) V4 - tan'^ \ - + log x 

V2 + X 



27ra;-3x2 + |7ra;3 + 0(a;^). (5) 



We use the Passarino-Veltman one-loop functions [16] defined as 
A(m) = -^^(m^)i-% 

e Jo 

1 /■! 

522(^»^ml,m2) = — -/ dz[{l- z)ml + zml- z{l- z)p'^ -i6Y~% (6) 

2e{l — ej JO 

and ^ 

522(/,mi,m2) = 522(p^ml,m2) - -[v4(mi) + A(m2)]. (7) 

The 0(my) term of Eq. ([5]) gives the nondecoupling mass splitting within the multi- 
plet. For example, for F = 0, the one-loop mass splitting between and the neutral 
component p of F is written as [13] , independent of J, 

M(/«)-M(/°) = g2AM«. (8) 



where, in the M ^ limit, 

AM(^) = y (mvK - c^m^) = (166.99 ± 0.07)MeV. (9) 

The numerical value in Eq. ([9]) is obtained by using the pole masses mw = (80.398±0.025) 
GeV, mz = 91.1876 GeV, 0:2 = a{'niz)/s^r = a;(m^)/(l— m^/m|), and the QED running 
coupling in the MS scheme a{mz) = (127.93 ± 0.03)~^, cited from Ref. [IT], as input 
parameters. Note that the value ([9]) should change by ~ 1 MeV depending on choices of 
the renormalization scheme for the input parameters. 



3 Two-loop mass correction 

We now calculate the two-loop mass correction SM^'^^ coming from diagrams shown in 
Fig. l(b-e). We use Feynman gauge fixing for simplicity, although the final result should 
not depend on the gauge fixing method. 

The contribution of the diagram Fig. 1(b) with the insertion of the one- loop self 
energy of the electroweak gauge boson, lV^l}''^{k) = g^^U^^^'^{k^) + 0{k^k^), is written as 

<5m(2'1) = -{CF-li)AEww-sl.{h + Y)^AJ:^^ 

-2swcw{h + Y)ih - tlrY)AE^z - c^wih - tl^YfAllzz, (10) 
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(a) 



(b) 





Or-, 



Figure 1: One- loop (a) and two-loop (b-e) contributions to the self energy of the heavy fermions 
/ in the multiplet F. The solid thick line and wavy line represent F and electroweak gauge bosons 
V = {'y,Z,W), respectively. The black circle in (b) represents the one-loop self energy of the 
gauge bosons II^i^^^ standard model particles. 



where 



AE 



V1V2 



{2n)D [P - m2,J[A;2 - m^. ][(A; + - M^] 



■P=M 



Here we list the analytic forms of 11^^ ^^(/c^) in the standard model for completeness [18] 
The contributions from the {t,b) quark loops are, up to the overall factor Nca2/{4:Tr 



{Nc = 3 is the color number 


of quarks), 






unk') = 


^ [-4B22{kl,mt, 


{])-{kl-ml)B^{klmt,Q) 


(12) 


^Vik') = 


q=t,b 


), 


(13) 








(14) 




q=t,b^W 




) 



n 2 



Bo{k'^,mt,mt) 



where 



(15) 
(16) 



U^e) = -8B22{klmq,mq)-2klBo{klmq,mq). 

The contributions of other quarks and leptons are obtained by appropriate changes of 
mt, Qq, and Nc- For the gauge and Higgs boson loops, we have, abbreviating the overall 



factor ail (47r), 



ill 



.2 



1 - e)522(A;^ mw, 0) + Ak^Boie, mw, 0)] 
+ (1 + 8(1 - e)4y)B22{k\ mw, mz) 
+ [4:Cwk^ + 3m^]So(/c^ mw, mz) 

+522(fc^ mw, mh) - m'^Bo{k^, mw, mn), 
s^[4(3 - 2e)B22{k'^, mw, m,w) + 4:k'^Bo{k'^, mw, mw)]. 



^2 



9„2 



B22{k^,mw,m,w) 



+swCwi4:k'^ + 2m|)i?o(fc^, mw, m,w) 
1 



(17) 
(18) 



(19) 



4c'w{3-2e)-A 



^w 



B22{k ,mw,mw) 



+ [4:C^k^ + 4m^ - 2m|]5o(P, mw, mw) 
1 ~ 

+—[B22{k^,mz,mh) - mlBo{k^ ,mz,mh)]. 

Cw 



(20) 



In addition, there are also the contributions of F to H^^^^. However, it is shown that the 
resulting 0{mw) contributions to SM^"^^ are completely cancelled by the renormalization 
of the parameters in SM^^^ . 

We may calculate the integrals ffTTl) by extending the general formulas for the two- 
loop mass corrections [19], by including finite masses for {W,Z). However, since we are 
interested in the SU(2)2,-breaking and nondecoupling part of Eq. flTTl) . it is prefered to 
expand the integrals (fTTj) in mw{^ mz,mt,mh) and then separate the 0{mw) terms 
from the dominant and gauge-symmetric 0{M) terms, before numerical evaluation. This 
is achieved by applying the asymptotic expansion of the Feynman integrals near the 
threshold = M^, as described in Ref. [20]. The 0{mw) part of the integral (ITTl) is then 
obtained as 



AS 



ViV2\0{mw) 



d^k 



2M 



(27r)« [P - mljk^ - ml^]{2k-p) 



H 



V1V2 



(P). (21) 



In the following, we show only the 0{mw) part fl^ of the corrections ASv^y^. By 
substituting the self energies (fT2] - [20l) . the integrals (!2T|) are expressed in terms of the 
two- loop functions (a = 1, 2) 



(47r)2 

i 



Xo^a{mv,mi,m2) 



(47r) 



X22-a{mv,mi,m2) = J 



d^k 




M 




{2ti)d 


[k^ 


— myY 


{2k -p) 


d^k 




M 




{2tx)d 


[k^ 


- ml^Y 


{2k -p) 



Bo{k'^,mi,m2) - 



B22{k^,mi,m2) 



1 I ml + ml 
~t \ 4 
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f22) 
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and products of the one-loop functions. Note that the functions in Eq. fl22|) are indepen- 
dent of M and have no overall divergences. We calculate these functions by numerical 
integration of the Feynman parameter integrals shown below, 



Xo_i(mv, mi, = vrmy 



log- 



m 



V 



dz 2 



ri(l - z) + r2Z 



z(l-z) 



- 2 log(Vri(l -z)+ r^z + J z{l - z)) 



Xo_2{mv,mi,m2) 



IT 



2mv 



log 



+ 2 f^dz \ 



(23) 



z(l-z) 



- z) + r2Z + Jz(l - z 



X22-i{mv, nil, 1712) 



+ log(^ri(l - z)+ r2Z + ^Jz{l - z)) 
"1 



(24) 



TT 

m' 



V 



9 



-{1 - 3(ri + ra)} log — - - -(ri + rs 



4 



+ dz { {-3z{l -z)+ 2(1 - z)ri + 2zr2} 



1 — 2;)ri -|- zr2 



z(l-z) 



+3{z{l - z) - {1 - z)ri - zr2} l^\og{^Jri{l - z) + r2Z + z{l - z)) 



--log2(l - z) 



(25) 



X22-2{mv,mi,m2) 



7T 

--my 



m 



1 - ri - rs) log — - - + -(ri + rg) 
/i"^ / 4 



+ d^; (^-3^z{l - z)^{l - z)ri + zr2 + {32:(1 - z) - (1 - z)ri - zrs} 



X |log(^(l - z)ri + 2;r2 + z{l - z)) - - \og{z{l - z)) 



1 
2 

where ri_2 = m\2/niy. fi is the MS renormalization scale. 

Here we show the explicit forms of the integrals (|2T|) . after subtracting 0(l/e) diver- 
gences from IlviV2 by the MS scheme, and separating the mass corrections to the gauge 
bosons 5my = —Kell^^ {■my){V = W, Z) from 11^^ [k"^). The {t,b) contributions coming 
from Eqs. ( |T2] - [T6l) are, up to the overall factor A^c«2/(47r), 



(26) 



Aytb 
2 Aytb 

2swCw^^^'^z 



TT 



2{tb) 
W ' 



1 

2 miv 

2slrQt{ht - 2Qtslr)X^z{mt) 
+2sl,Qb{hb - 2Qhsl,)X^z{0), 



(27) 
(28) 

(29) 
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c^AS*J^ = [{htY - 2htQtsiv + 2Qls%,\Xzz{mt) - mfGoimz, rrit, rrit) 

+ [{hb)^ - 2hbQbslr + 2QlsUXzz{0) - cl.—6mT\ (30) 

mz 

with 

Xwy/irrit) = 8G22{mw,mt,0) + 2Xo.i{mw,mt,0) 

+2{mly - m^)Go{mw,mt,0), (31) 

X^zimt) = -^X22-i{mz,mg,mg) ^ 2" 

+2Xo^i{mz,mt,mt) + 1 h _ log ^ , (32) 

mz \ J 

/4, ml 20\ , , 

X,z{Q) = 7rm^Ulog-f --j, (33) 

Xzzimt) = 8G22imz,mt,mt) 

+2Xo^i{mz,mt,mt) + 2mlGo{mz,mt,mt) (34) 

XzziO) = ^mzl^Uog'^-^y (35) 

Here we used the notations 

G22{mv-,mi,m2) = X22-2i'mv,mi,m2) 

+ - [Re522 (my, mi, m2) - ((m^ + m2)/4 - my/12)/e], 

Zmy 

vr 

Go{mv,mi,m2) = Xo_2{mv,mi,m2) + [Re5o(my, mi, m2) - 1/e], (36) 

Zmy 

and substituted analytic forms of the two-loop integrals (122|) at mi = m2 = and at my — > 
0. Analytic forms of other integrals involving mt are shown in Appendix. Contributions 
of other quarks and leptons are obtained by taking m^ — > 0, where 



XH.w'(0) = vrmH^(^log^-^), (37) 



and, for leptons, changing {Qg, N^). Similarly, the gauge and Higgs boson contributions 
coming from Eqs. (fT71 - [20l) are, up to the factor 0:2/(47?), 



AEJ^V = -16s^G22{mw,mw,0) - 8s'^{Xo_i{mw,mw,0) + mlyGo{m'^,mw,0)) 

r + l)G22(miy,mvi/,?T 

-2((5m^ + m?z)(?w - ^zSw)Go{mw, m^y, mz] 



-2{8cl^r + l)G22{mw, mw, mz) + -nmw - 8c^Xo_i(mvp', m^y, mz) 



-2G22{mw, mw, mt) + 2m'^Go{mw, mw, mu) —5m^^^\ (38) 

mw 

sI^^TX^!; = lOTT^stymiy, (39) 
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Vh 



mr 



-X22-i{mz,mw,mw) 



2 2 S 



IT 



m 



AS 



4 

-2(12c^ - 4c^ + l)G'22("^z, "^ly, mw) + ■^c^'^mz 
-8c^Xo_i{mz,mw,mw) - 4:{4:C^ - l)m^[Go{mz,mw,mw)] 



(40) 



TT 



-2G22("^z, mz, rrih) + 2m^Go("^z, "^z, "^/i) - c,^ 6m 



2{Vh) 



mz 



(41) 



Other diagrams shown in Fig. l(c-e) are also evaluated by using the threshold expan- 
sion [20], keeping only the 0(my) parts. Their sum, with subtracting sub divergences by 
the MS scheme and after the (one-loop) x (one-loop) term in Eq. ([3]) is added, is given as 

i2 



c"^log^ + (2-c^^)log 



m 



w 



-l + log2) 



m 



+87rciymzhih - t'wY) log - ^T^4,{Cp - Ii)fzw- 



(42) 



Here 



zw 



1 

"3' 
1 

3' 

-0.027m 



;(2 + c^w)mw 
-^(2 + c^)mz 



Jo 

dzz-^/\l-. 



-1/2 



[{c^'z + l-zfl^-l] 



1-1/2 



[{cl^z + 1 



N3/2 



1] 



(43) 

is the two-loop function appearing in Fig. l(c,d) with both W and Z bosons. 

We then need to add the counterterms coming from the renormalization of the pa- 
rameters in the one-loop contributions (jH [5]); (mvi/,m^) in Xy^^z (0^2, c^, ...) in 
the coupling constants. We adopt the scheme where the pole masses {mw, mz) and the 
MS running coupling of QED a(m^), which are used in Eq. ([9]), are chosen as the input 
parameters. In this scheme, the renormalization is achieved by removing the last 0{5m\r) 
terms from ASvi/vi^ (pTl [38l) and AE^^ (j30ll4TI) . and adding the counterterms for (0^2, c^, 
. . . ) expressed as tree-level functions of (mz, mw, Oi{mz)). It is checked that the final 
form of the two- loop 0{mw) mass correction to is finite and independent of the MS 
renormalization scale /i. 

Here we comment on the mass splitting of a new heavy scalar SU(2)l multiplet S. In 
contrast to the case of the fermion multiplet, direct couplings of S to the Higgs bosons, 
such as S*SH^H, should always exist [T3l[2T]. Nevertheless, assuming that the effect of 
these direct couplings is negligible, we have verified that the nondecoupling 0{mw) parts 
of the one-loop 122] and two-loop mass corrections 5M are identical to those for the 
fermions in the same gauge representation. This result is quite natural in the view that 
the 0{mw) mass correction could be understood as the energy of the electroweak gauge 
fields around a static point source, and should be insensitive to the spin of the source 
particle [T3] . 
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4 Numerical results 



We show the numerical results of the two-loop contributions to the mass splitting within 
the Y = fermion multiplet. As seen in Eqs. fllO[ the one-loop relation ([H]) still holds 
with the change AM^^) ^ AM^^) + AM^^\ where AM^^) = AM^^-sO + AM^^-^^). 

The contribution AM'^^''^') of the quark-lepton subloop diagrams (including corre- 
sponding counterterms) is shown in Fig. 2 as a function of m^. At rrit = 171 GeV, 
there is cancellation between the {t, b) subloop contribution, shown in the dashed line, 
and remaining contribution with subloops of other quarks or leptons, —3.3 MeV, giving 
the total shift -1.5 MeV at mi = 171 GeV. 




.3 I I I I I I I 

140 150 160 170 180 190 200 

mt [GeV] 



Figure 2: Two-loop contribution to the mass splitting AAf^^''^^^ between fermions in a heavy 
SU(2)i multiplet with y = 0, from diagrams in Fig. 1(b) with quark and lepton subloops. Solid 
and dashed lines denote full and {t, h) subloop contributions, respectively. 

The remaining contribution AM^^'^'^^ from diagrams without quarks or leptons (again 
including corresponding counterterms) is shown in Fig. 3 as a function of 171^. At nth = 
140 GeV, the shift is —0.9 MeV, smaller than the quark-lepton loops. 

These two- loop contributions are much smaller than the 0{mw) part of the leading 
one-loop contribution (jl]), as expected. However, for F = 0, it may compete with the 
M-dependence of the one-loop contribution ([8]) which behaves like —0.5(1 TeV/M)^ MeV 
for large M due to the accidental cancellation of the 0(m^/M) term in Eq. ([8]). In 
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Figure 3: Two-loop contribution to the mass splitting AM^'^'^^^ between fermions in a heavy 
SU{2)l multiplet with Y = 0, from diagrams in Fig. l(b-e) with gauge and Higgs bosons. 

comparison, in the case of the higgsino-like doublet F = (/^,/") with (/ = 1/2, Y — 
1/2), the two-loop corrections to the mass splitting M{f^) — M(f^), which is amz/2 = 
356.4 MeV at the one-loop, is —1.2 MeV from quark and lepton loops at rrit = 171 GeV, 
and —1.8 MeV from gauge and Higgs boson loops at rrih = 140 GeV, respectively. 

5 Conclusion 

We have calculated the two-loop electroweak contribution to the 0('mw) correction to the 
masses of new heavy fermions in an SI]{2)l multiplet F, which causes gauge-symmetry- 
breaking and "nondecoupling" mass splitting within F. Analytic formula of the 0{mw) 
mass corrections have been presented for F in general 811(2)^x11(1) representation. The 
two- loop contribution has turned out to be typically O(MeV), which is of similar order 
to the M dependence of the one-loop contribution for the Y — case. 
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Appendix 



In the case of the {t, h) contributions flTTHSUI) . Feynman parameter integrals for the func- 
tions (I22p can be analytically performed. For mi = m2 = ^frrny with r > 1/4, we 
have 



-^0-1 

X22-1 



TTinv 
71 



log 



2 + log r - 2 V4r - 1 tan ^ ^Ar - 1 



log- 



2m\/ 

vr 

n 

3 

1 — 6r 



+ log r + 



V 



m 



V4r - 1 
^ 2 7 



tan ^ V4r - 1 



(A.l) 
(A.2) 



-(1 - 6r) log ^ o + o'^ 



3 2 



^ logr + -(4r - l)^/^tan ^ V4r - 1 



22-2 



TT 



m 



^(l-2r)log^-^ + ^r 



1 



4 



2 2 



1 - 2r , 1 , , 

H — logr v4r — Itan v4r — 1 

4 2 



(A.3) 



(A.4) 



For 7712 = and mi = y^my, 



-^0-1 
-^0-2 
-^22-1 



vrmy 

TT 

2my 

^ 3 

mV 



log 



m 



V 



2 — 2^/r + r log r — 2(r — 1) log(l + 



771 



log — ^ — 2v^ — r logr + 2(r + 1) log(l + ^/r) 



12- 



(A.5) 
(A.6) 



-(1 - 3r) log ^ + — (-8 - 6v^ + 21r + I6r^/^ - + Gr^/^) 
4 ^ 1 2 



('^ _ 3y2 ^ 
+ - logr - -(r - 1)=^ log(l + v^) 



(A.7) 



X. 



TT 



22-2 



-my 



2 

r^(r — 1) 



-(1 - r) log ^ + —(-6 - 6Vr + 3r + 4r^/^ + 3r2 - Gr^/^) 
4 /i^ 12 



^ logr + -(r - l)(r^ - 1) log(l + v^) 



(A.8) 
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